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Protein function often involves changes between different conformations. Central questions are how these 
conformational changes are coupled to the binding or catalytic processes during which they occur, and how they 
affect the catalytic rates of enzymes. An important model system is the enzyme dihydrofolate reductase (DHFR) 
from E. coli, which exhibits characteristic conformational changes of the active-site loop during the catalytic 
step and during unbinding of the product. In this article, we present a general kinetic framework that can be used 
(1) to identify the ordering of events in the coupling of conformational changes, binding and catalysis and (2) to 
determine the rates of the substeps of coupled processes from a combined analysis of NMR R2 relaxation disper- 
sion experiments and traditional enzyme kinetics measurements. We apply this framework to E. coli DHFR and 
find that the conformational change during product unbinding follows a conformational-selection mechanism, 
i.e. the conformational change occurs predominantly prior to unbinding. The conformational change during the 
catalytic step, in contrast, is an induced change, i.e. the change occurs after the chemical reaction. We propose 
that the reason for these conformational changes, which are absent in human and other vertebrate DHFRs, is 
robustness of the catalytic rate against large pH variations and changes to substrate/product concentrations in 
E. coli. 



INTRODUCTION 



The structural changes in an enzyme are intimately linked 
to its function. Ligand-exchange events of enzymes through 
either binding/unbinding or chemical transformation are of- 
ten accompanied by conformational changes to the protein 
0]|2]. Indeed, the rates of catalysis are frequently limited 
by the slow conformational changes in enzymes (e.g. SQ). 
Optimal enzyme engineering would thus require a deep un- 
derstanding of both the chemical steps of catalysis and how 
the structural changes in the enzyme help to facilitate these 
events j6]-[8). Of course, one powerful source of functional 
information comes from the study of enzyme kinetics. How- 
ever, classical enzyme kinetics have traditionally followed the 
chemical transformation of substrate to product, and only give 
indirect evidence about structural changes to the enzyme it- 
self. For example, substrate binding involves the formation of 
noncovalent interactions between enzyme and substrate, but a 
conformational change in the enzyme (and/or substrate) may 
be required to maximize these interactions; traditional enzyme 
kinetics may not be able to tease out these two separate, but 
related events. More recent biophysical techniques, including 
those from single-molecule spectroscopy 191 4TTI and solution- 
state nuclear magnetic resonance (NMR) |[T2HT6l . can moni- 
tor structural changes in an enzyme itself more directly. A 
synthesis of kinetic information from these recent biophysical 
experiments and from more traditional ensemble enzyme ki- 
netics would bring a richer understanding of the physical and 
chemical barriers to catalysis. 

One technique that has brought deep insights into the "con- 
formational gymnastics" an enzyme undergoes during its cat- 
alytic cycle is NMR R2 relaxation dispersion spectroscopy 
|[T2l IT3l [161 [171 . This technique measures conformational ex- 
change between two (or more) protein conformations, and can 
give insight into the kinetics of exchange, the populations of 
each of the exchanging conformations (i.e. the thermodynam- 
ics), and can give structural information about higher energy 



conformations, even when the conformational equilibrium is 
highly skewed (i.e. less than 5% of the conformational en- 
semble is in the higher energy conformation) |[T4l [161 [17). 
Such methodology has been used to monitor conformational 
exchange processes on the /is-ms timescale for all of the inter- 
mediates in the catalytic cycle of Echerichia coli dihydrofo- 
late reductase (DHFR) 0(18). DHFR catalyzes the reduction 
of dihydrofolate (DHF) to tetrahydrofolate (THF) using the 
cofactor NADPH, and has served as an important model sys- 
tem for investigating the interplay of conformational dynam- 
ics, binding and catalysis |[T9]j26l . Intriguingly, the higher 
energy or "excited-state" conformations for all the intermedi- 
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FIG. 1: E. coli DHFR cycles between closed (left; PDB 1RX2) and 
occluded (right; PDB 1RX4) conformations. These conformations 
are dependent on interactions between three loops. In the closed con- 
formation, the active-site Met20 loop (in blue) hydrogen bonds to the 
FG loop (in red). Attainment of the occluded conformation requires 
breaking of these hydrogen bonds and formation of new hydrogen 
bonds between the Met20 and GH (in purple) loops. 
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ates in DHFR catalysis are structurally similar to either the 
next and/or previous "ground-state" conformation in the cat- 
alytic cycle of DHFR [4|. In other words, the excited-state 
conformations predict the next structural state of the enzyme. 
With such rich functional information, E.coli DHFR is ideal 
in the development of a kinetic framework that encompasses 
both the traditional pre-steady state enzyme kinetics and the 
kinetics on the structural changes in the enzyme gleaned from 
the NMR i?2 relaxation dispersion experiments. 

The i?2 relaxation dispersion results from DHFR also 
prompted the suggestion that ligand-exchange events on this 
enzyme, and likely other enzymes lfl6l . can operate through 
a "conformational selection" mechanism. In conformational 
selection l27l . conformational changes occur predominantly 
prior to a binding process, an unbinding process, or a catalytic 
process; the ligands appear to select a conformation from the 
protein ensemble for binding, unbinding, or catalysis fl4ll28T 
[36]. In contrast, "induced change" |37| proposes that the con- 
formational change occurs predominantly after a binding, un- 
binding, or catalytic process [38-41]. Critically, the timing 
of the ligand interaction and the structural changes to the en- 
zyme is not trivial; the processes of conformational selection 
and induced change have important biological consequences 
in terms of substrate promiscuity, enzyme regulation, protein 
evolution, and biological information storage Bl l42l l43l . 

In this article, we present a general kinetic framework that 
can be used (1) to distinguish conformational-selection from 
induced-change processes, and (2) to determine the rates of 
the physical substeps (conformational changes) and chemical 
substeps (binding/unbinding, or chemical reactions) of these 
processes from a combined analysis of NMR i?2 relaxation 
dispersion and traditional enzyme kinetics experiments. We 
apply this framework here to E. coli DHFR and derive ex- 
tended catalytic cycles of this enzyme that specify the or- 
dering of events in the coupling of conformational changes, 
binding and catalysis. We focus on the major conformational 
changes of E. coli DHFR that occur during the catalytic step 
(from "closed" to "occluded", see fig. [TJ and during unbinding 
of the product THF (from "occluded" back to "closed"). We 
find that product unbinding occurs via conformational selec- 
tion, whereas the conformational change during the catalytic 
step is an induced-change process (see Results section). Our 
extended cycles help to explain the role of the conformational 
changes for the function of E. coli DHFR (see Discussions 
section). 

The methods employed here to analyze the catalytic step 
and product unbinding of E. coli DHFR can be generalized 
to binding, unbinding, and catalytic processes of other protein 
systems. Our methods are based on general expressions for 
the effective rates of processes that involve two or more sub- 
steps, and on a mutational analysis of these effective rates. In 
this analysis, we focus on mutations distal to the binding site 
that mainly affect the conformational equilibrium of the pro- 
tein. We find that conformational-selection processes are sen- 
sitive to such distal mutations because these processes involve 
a change to a low-populated, excited-state conformation, and 
because the equilibrium probability and excitation rate of this 
conformation depend on the conformational free-energy dif- 
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FIG. 2: Catalytic cycle of the enzyme dihydrofolate reductase 
(DHFR) fromii. coli (48 1. DHFR catalyzes the reduction of dihydro- 
folate (DHF) to tetrahydrofolate (THF), using NADPH (NH) as a co- 
factor. The enzyme (E) cycles through five major intermediates! 19. 
I20ll48l . X-ray crystal structures of the enzyme in these five inter- 
mediates indicate characteristic conformational changes [45 1. In the 
intermediate states E NH and Eqhf, the active-site loop is predomi- 
nantly 'closed' over the reactants (see fig[TJ. In the other three in- 
termediates, this loop 'occludes' (i.e. protrudes into) the active site. 
The conformational changes occur during the catalytic step, and dur- 
ing unbinding of the product THF. In this so-called pH-independent 
scheme [48], the forward rate of the catalytic step is the maximal rate 
at low pH, and the reverse rate of this step is the maximal rate at high 
pH (see text). 



ferences. In contrast, induced-change processes involve a con- 
formational relaxation into a new ground-state after a binding, 
unbinding or catalytic event, which is rather insensitive to 
changes in conformational free-energy differences, provided 
(i) the conformational relaxation is fast compared to the rates 
of the chemical substep PBl . or (ii) the transition-state for the 
conformational exchange is close in free energy and structure 
to the excited protein conformation. The analysis of the effect 
of distal mutations thus can provide the basis for a simple di- 
agnostic to identify conformational-selection versus induced- 
change processes. 



RESULTS 

Structural changes in the catalytic cycle of E. coli DHFR 

The catalytic cycle of E. coli DHFR involves five major 
intermediates (see fig.|2]i. X-ray crystal structures indicate an 
'occluded' ground-state conformation of the active-site Met20 
loop (residues 9-24) in the three intermediates with bound 
product THF, and a 'closed' conformation of this loop in the 
other two intermediates B31 . When the Met20 loop is in an 
occluded conformation, its central region (residues Metl6 and 
Glul7) is flipped into the cofactor binding pocket such that the 
nicotinamide ring of the cofactor is sterically occluded from 
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FIG. 3: Extended catalytic cycle of the enzyme DHFR from E. coli. Excited-state conformations are shown in red, ground-state conformations 
in blue. On this extended cycle, the catalytic step is decomposed into the actual chemical substep in the closed conformation of the enzyme 
required for hydride transfer, and a physical substep in which the enzyme conformation changes from closed to occluded. The given rates 
for the chemical substep follow from eqs. |4]l and {5} and the experimentally determined maximal rates for the overall catalytic step shown in 
fig. [2] The transition rates between the ground-state and excited-state conformations have been determined in NMR relaxation experiments 
t4l . Our analysis of the product-unbinding kinetics and NMR relaxation experiments |4| indicate that the conformational change from the 
occluded to the closed conformation occurs largely prior to the unbinding of the product THE Along our extended cycle, the product unbinds 
from an excited state with a conformation similar to the closed conformation (denoted as 'closed^'). The given binding and unbinding rates 
of THF follow from eqs. l |19} and \2Q) and the overall rates for the product binding and unbinding process shown in fig. [2] 



binding. However upon binding cofactor, residues Met 16 and 
Glul7 are turned away from the nicotinamide binding pocket, 
residues Metl6-Alal9 zip up into a type III' /3-hairpin, and the 
side-chains of Asnl8 and Met20 close down over the pterin 
ring of the substrate DHF to form the closed conformation. 
Only the closed conformation allows proximal positioning of 
the cofactor and substrate reactive centers and therefore it is 
the closed conformation that must be adopted in the Michaelis 
complex. The closed and occluded Met20 loop conformations 
are stabilized by unique hydrogen bonding networks to the FG 
(116-132) and GH (142-150) loops, respectively (see fig. [T). 
Changes to the amino acids in the FG or GH loops that make 
contact with the Met20 loop can severely alter the equilibrium 
between the closed and occluded conformations of the Met20 
loop. The mutant G121 V, for example, has the occluded con- 
formation as ground state in all five intermediates along the 
catalytic cycle [46], while the mutant S148A seems to have 
the closed conformation of the active-site loop as ground state 
in all intermediates (47). These amino acid substitutions have 
functional consequences for both the chemical step and the 
THF product release step that are accompanied by closed- 
occluded conformational changes. 



The catalytic step of E. coli DHFR operates through an 
induced-change mechanism 

Two steps in the catalytic cycle of E. coli DHFR are ac- 
companied by major structural changes in the enzyme - the 
catalytic step itself and the release of the product THF from 



the ternary complex of the enzyme with THF and NADPH. 
Division of the catalytic cycle into chemical (including sub- 
strate binding/product release and chemical transformation of 
substrate to product) and physical substeps requires an inves- 
tigation of both steps. Hydride transfer in the chemical step is 
only possible in the closed conformation, thus any conforma- 
tional change in the enzyme through the catalytic step must 
occur after THF production i.e. through an induced-change 
process. In the following, we apply our kinetic framework to 
determine the rates of the chemical and physical substeps of 
this induced-change process from experimental data. 

The catalytic step of wild-type DHFR is associated with a 
conformational change from the closed conformation (cE) to 
the occluded conformation (oE) of the enzyme (see fig. |2j: 

cE DHF Y ° E THF W 

The forward and reverse rates kf and k r for the catalytic step 
are pH-dependent since the catalyzed reaction DHF + NH + 
H + ?± THF + N + requires an additional proton H + in the 
forward direction. The forward rates k / from enzyme kinetics 
measurements at different pH are well described by ||48| 

k f (pH)^k md J(l + 10" H -^) (2) 

with the maximal forward rate fc max — 950 ± 50 s _1 at low 
pH and pKa = 6.5 ± 0.1 at 25°C. The reverse rate k r attains 
its maximal value 0.55 ±0.1 s _1 at high pH [48 1 and has the 
value k r ^ 0.03 s" 1 at pH 7 l49l . 

We now break down the catalytic step into its chemical and 
physical substeps. Since the hydride transfer from NH to DHF 
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is only possible in the closed conformation (see above), the 
chemical hydride-transfer substep occurs before the physical 
substep, the conformational change from the closed to the oc- 
cluded conformation: 



cE^ 



k+ 



cE* 



r ' or 



(3) 



The occluded conformation oE^f is the ground state of the 
enzyme with bound products THE and N + , and the closed 
conformation cE^ F is an excited state. The rates k oc and 
k co for this conformational transition have been measured in 
NMR i?2 relaxation dispersion experiments. At 27°C, these 



rates are k n 



18.5 ± 1.5 s- 1 and k c 



1280 ± 50 s" 



The i? 2 relaxation dispersion experiments indicate that 
the conformational dynamics of the active-site loop does not 
depend on the pH l50l . 

The rates k + and fc_ for the chemical substep in the closed 
conformation can be calculated from (i) the rates k r and kf 
for the overall catalytic step determined in enzyme kinetics 
experiments, and (ii) the rates k oc and k co for the physical 
substep determined in i?2 relaxation dispersion experiments. 
Our calculations are based on general results for the effec- 
tive rates of processes with two (or more) substeps (see Ap- 
pendix). We first focus on the reverse direction. The effective 
rate k r for the reverse hydride transfer depends on k oc , k co 
and k- . Since the excitation rate k oc is much smaller than the 
ground-state relaxation rate k co , the effective reverse rate is 
k r ~ fc oc fc_/(fc co + (see eq. (32i in Appendix), which 
leads to 



-kr ^ (69 ±6) • k r 



(4) 



Product release from E. coli DHFR may occur through three 
potential routes 

Product release is the rate-limiting step in the steady-state 
catalytic cycle of DHFR for pH values up to the preferred 
internal pH of E. coli (see Discussion section below), and 
so, has a special importance in DHFR function. In contrast 
to the catalytic step, the unbinding of the product may oc- 
cur along three possible routes on which the conformational 
change from the occluded to closed conformations occurs ei- 
ther (1) prior to, (2) after, or (3) partly prior to and partly af- 
ter the unbinding event (see Figures 4 and 5). By calculating 
the effective rates along these three routes and by determining 
the effect of distal mutations on the off -rates, we can identify 
which of these routes is the product unbinding route of E. coli 
DHFR. It should be kept in mind that these three routes can 
be generalized to other protein systems and are applicable to 
any ligand (small molecule and/or large macromolecule) bind- 
ing/unbinding event, including binding events that involve 
folding of a ligand ||5Tll52l . An implicit assumption of these 
routes is that the chemical and physical substeps can be sepa- 
rated in time. 

Route 1: We first consider the conformational-selection 
unbinding route 
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c+[THF] 



CE^ 



(6) 



along which the conformational change precedes product un- 
binding. We assume here and below that the concentration 
of the product THF is significantly larger than the concentra- 
tion of the enzyme, which implies pseudo-first-order kinetics 
with approximately constant product concentration and, thus, 
constant binding rate c + [THF] . 



since k oc is much larger than k r . We have inserted here 
the values for the conformational transition rates k oc and 
k co measured in the NMR R 2 relaxation dispersion experi- 
ments. For the experimentally determined maximal effective 
rate k r = 0.55 ± 0.1 s _1 at high pH, we obtain = 38 ± 8 
s _1 as maximal reverse rate. 

The effective rate kf in the forward direction depends on 
the rates fc + , and the relaxation rate k co into the ground 
state with bound products THF and N + . Because k co is much 
larger than fc_, the effective rate in the forward direction 
is either k + or k co , depending on which of these two rates 
is smaller and, thus, rate-limiting (see eq. (32l). Since the 
ground-state relaxation rate k co is larger than the maximum 
effective rate in the forward direction, we have 



route 1 : conformational selection 



closed 



4 




Q+[THF] w /£nh\ 
^ ^TH^ 



occluded 



(5) 



This equation is consistent with eq. Q because the equilib- 
rium condition kf/k r ~ (fc oc //c co )(fc+/fc_) is fulfilled. Our 
decomposition of the catalytic step into its physical and chem- 
ical substeps is included in the extended catalytic cycle of 
E. coli shown in fig. [3] 



route 2: induced change 

FIG. 4: Possible unbinding routes of the product THF from wildtype 
E. coli DHFR. On the conformational-selection route 1, the confor- 
mational change from the occluded to the closed conformation of the 
active-site loop occurs prior to product unbinding. On the induced- 
change route 2, the conformational change occurs after unbinding. 
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route 3: unbinding from intermediate conformation 




closed intermediate intermediate occluded 



FIG. 5: Possible unbinding route of the product THF from wildtype E. coli DHFR. Along this route, the product unbinds from an intermediate 
conformation. 



The effective off-rate fc^jj along route 1 depends on the 
rates b oc and b co of the conformational transitions and on the 
unbinding rate c_ from the closed conformation. Since the 
excitation rate b oc is much smaller than the relaxation rate b co 
into the product-bound ground state, the effective off-rate is 



(see eq. (32 1) 



1.(1) 



b co + c_ 



(7) 



The effective on-rate constant kon on route 1 depends on the 
binding rate constant c + and the unbinding rate c_ in the 
closed conformation, and on the relaxation rate b co into the 
product-bound ground state. If the relaxation rate b co is much 
larger than the binding rate c + [THF] , the effective on-rate con- 
stant is (see eq. d32ji) 



fen* 



b co + C- 



(8) 



The eqs. (jTjl and <|8j fulfill the equilibrium condition 

k$/k$ = (6 co /6 oc )(c+/c_). 

Amino acid substitutions in the FG and GH loops of E. coli 
DHFR, such as G121V and S148A, affect the conformational 
equilibrium between closed and occluded states, but these 
residues do not make direct interactions with the THF prod- 
uct or NADPH cofactor themselves. These distal mutations 
thus should mainly affect the free-energy difference AG co be- 
tween the closed and occluded conformation. We will show 
in the following that the analysis of distal mutations in the 
context of the three routes of product unbinding can be used 
to separate out the chemical and physical substeps of product 
release and give insight into the ordering of the events. 

In transition-state theory, the conformational transition 
rates b oc and b co depend on free-energy differences to the tran- 
sition state: 



6 -AGW«r and b 



b Q e 



-AG„ a /RT 



(9) 



Here, AG* is the free energy difference between the transi- 
tion state and the occluded conformation, and AG* C the free- 
energy difference between the transition state and the closed 
conformation. Because of b oc -C b co , it seems not unrea- 
sonable to assume that the free-energy difference AG* C is 
significantly smaller than AG* , i.e. that the transition state 
is significantly closer in free energy to the closed conforma- 
tion (see fig. [6]). According to the Hammond-Leffler postu- 
late, the structure of the transition state then resembles the 



structure of the closed conformation |f53~ll54l . Distal mutations 
thus should have a similar effect on the transition state and the 
closed conformation, which implies AAG* G ~ AAG co and 
AAG* C <C AAG co for the mutation-induced changes of the 
free-energy differences. This leads to 



b'oc/bo 



-AAGco/ RT and y 



(10) 



where the prime indicating rates of the mutant, and A AG co = 
AG^, — AG co is the mutation-induced change of free-energy 
difference between the conformations. We have assumed here 
that the pre-exponential factor b a in eq. |9]) its not affected 
by the mutation. According to eq. ( fT0| ), mutational shifts of 
the conformational free-energy difference affect the excitation 
rate b oc , but not the relaxation rate b co . 

If we neglect the effect of the distal mutations on the un- 
binding rate c_ from the closed conformation, we obtain 



i off 
"■off 



^— AAGco/ RT 



(ID 



from eqs. (|7]i and (lOi. In the special case c_ <C b co of a 
conformational relaxation rate b co that is much larger than the 
unbinding rate c_, eq. ( TP!) also follows directly from eq. (j7]) 

and b oc /b co = e~ AGco/ " for c'_ ~ e_, because we have 

Ktff ~ {b oc /b co )c_ in this case Pffl . However, this special 
case does not seem to apply to DHFR product unbinding (see 
below). 



Eq. (11 



simply states that the ratio of the mutant 
and wildtype off-rates for the conformational-selection pro- 
cess depends on the shift A AG co of the free energy-difference 
between the closed and occluded conformations induced by a 
distal mutation. 

Route 2: On the induced-change unbinding route, the con- 
formational change occurs after the unbinding of the product 
THF: 



oE NH 

THF 



o + [THF] 



oE 



NH 



cE 



NH 



(12) 



Along this unbinding route, the product THF unbinds from 
the occluded ground-state conformation, which induces the 
conformational change into the closed conformation. 

The effective off-rate rate for the induced-change route 
2 depends on the unbinding rate o_ from the occluded con- 
formation, the rebinding rate o + [THF], and the relaxation rate 
u oc into the unbound ground state. If the relaxation rate u oc 
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is much larger than the rebinding rate o+[THF], the effec- 
ts 

tive unbinding rate fc^ ff is either o_ or u oc , depending on 
which of the two rates is smaller (see eq. ([32])). Because 
of u co -C u oc , we assume again that the transition state for the 
conformational exchange is significantly closer in free energy 
and structure to the excited-state conformation (see fig. [6}, 
which implies that distal mutations should mainly affect the 
excitation rate u co rather than the relaxation rate u oc , as in 
eq. ( jlo| . Since the effective off-rate k^g does not depend on 
the excitation rate u co , we obtain 



k' {2) 

ft off 

J?) 

ft off 



(13) 



if we neglect the effect of the distal mutation on the binding 
rate constant o + and unbinding rate o_ in the occluded con- 
formation. 

In contrast to the conformational-selection process, the lig- 
and off-rate for the induced change-process is not affected 
by a distal mutation, according to eq. (13 1. The eqs. (Hi 
and ( 13 i thus provide the basis for a simple diagnostic to 



identify conformational-selection versus induced-change pro- 
cesses. In the next section, we will apply this diagnostic to 
experimental data for the product unbinding rates of wildtype 
and mutants of E. coli DHFR from the ternary complex E^ F , 
in combination with calculations of mutation-induced stabil- 
ity changes (see fig. 17}. We focus here on off-rates because 
corresponding data for on-rates are not available for DHFR. 
However, a corresponding diagnostic also holds for mutation- 
induced changes of on-rates ll44ll . 

Route 3: Our diagnostic to identify conformational- 
selection and induced-change processes can be extended to 
route 3 on which the product unbinds from an intermediate 
conformation i between the occluded and the closed confor- 
mation: 



b i 



oE? 



iE- 



NH 
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NH 



cE^ 



(14) 



ho 



i+[THF] 



Along this route, a conformation-selection step from the 
bound ground state oE^ F to the intermediate state iE NH is 
followed by an induced-change step into the unbound ground 
state cE NH . The routes 1 and 2 can be seen as limiting cases of 
route 3 with 'intermediates' that are identical with the closed 
or occluded conformation, respectively (see also |55 , 56 1). 



The effective unbinding rate of route 3 is (see eq. ( 39 1): 



1.(3) 
^off 



bnd- 



b io + i_ 



- i+[THF]) + u.; c z_ 
for u lc > i+[THF] 



(15) 
(16) 



For a relaxation rate Ui C into the unbound ground state that is 
much larger than the rebinding rate i + [THF] of the unbound 
excited state, the unbinding process thus is dominated by the 
conformational-selection step since fc^ ff is identical with the 
effective rate from the bound ground state oEjp/ F to the inter- 
mediate state iE NH in this case (see eq. for comparison). In 
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FIG. 6: Energy landscapes for the conformational changes along 
route 1 and route 2 of fig.|4]in transition-state theory. The full lines 
indicate free-energy levels for the wildtype, the dashed lines the free- 
energy levels for a mutant. Along route 1 , the conformational change 
occurs in the product-bound state in which the occluded conforma- 
tion is lower in free energy than the closed conformation. In the land- 
scape for route 1, the free-energy difference between the transition 
state and the closed conformation is small compared to the difference 
between the transition state and the occluded conformation. Accord- 
ing to the Hammond-Leffler postulate, the structure of the transition 
state then resembles the structure of the closed state, which implies 
that mutations have a similar effect on these two states. Along route 
2, the conformational change occurs in the unbound state in which 
the closed conformation is lower in free energy. Since the transi- 
tion state now is much closer in free energy to the occluded state, 
mutations have a similar effect on these two states according to the 
Hammond-Leffler postulate. 



analogy to eq. ( 1 1 1 for conformation-selection unbinding, the 



effect of a distal mutation on the unbinding rate therefore can 
be characterized by 



1.(3) 
""off 



-AAGia/RT 



for 



>?:+[THF] (17) 



where AAG io is the mutation-induced change of the free- 
energy difference between the intermediate and the occluded 
conformation. On route 3, the ratio of the mutant and wildtype 
off-rates thus depends on the shift AAG; of the free-energy 
difference between the intermediate and the occluded confor- 
mation. 
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TABLE I: Calculated mutation-induced stability changes for DHFR conformations 



state 



Ethf 
Ethf 



PDB file 


active-site loop 


D122A D122N D122S 


G121V M42W S148A 


1RX1 


closed 


0.22 


-0.08 


-0.52 


1.13 


-1.58 


0.34 


1RX2 


closed 


-0.28 


-0.79 


-0.45 


2.15 


0.41 


-0.26 


1RX4 


occluded 


-0.42 


-0.09 


-0.25 


-0.26 


-0.01 


1.04 


1RX5 


occluded 


-1.15 


-0.63 


-1.22 


-0.74 


-0.20 


0.31 


1RX6 


occluded 


-1.03 


-0.80 


-0.79 


-0.85 


-0.16 


0.79 



Stability changes A AG for DHFR mutants in units of kcal/mol calculated with the program Concoord/PBSA |60|. The calculation error can 
be estimated as the standard deviation of 1.04 kcal/mol between experimentally determined and calculated stability changes for 582 
mutations of seven proteins |60|. The given values are averages over five independent calculations with Concoord/PBSA. The standard 
deviation for the five calculations is on average 0.38 kcal/mol. The mutation-induced stability change is the difference AAG = AG' — AG 
between the stability AG' of the mutant and the stability AG of the wildtype. The stability AG is the free energy difference between a 
folded, native conformation and the denatured state. Atoms of the ligand molecules are discarded from the PDB files prior to the calculations. 



Mutational analysis of DHFR indicates that THF release 
operates through a conformational-selection mechanism 

Based on the above analysis, we can identify whether a 
ligand-exchange event operates through a conformational- 
selection or induced-change process by analyzing experimen- 
tal data for the ligand off-rates (or on-rates) for a series of 
enzyme mutants, in which the amino acid substitutions pri- 
marily alter the free-energy difference between the conforma- 
tions before and after binding. Benkovic and co-workers have 
measured the effect of six single-site mutations distal to the 
binding site on the product unbinding rate of DHFR ll49l [571 — 
59). Four of the mutations decrease the product unbinding 
rate, while two of the mutations increase the unbinding rate 
(see table [II]). Since most of the mutated residues are in con- 
tact with the active-site loop, the effect of the distal mutations 
on the unbinding kinetics seems to result predominantly from 
a shift of the conformational equilibrium between the closed 
and occluded conformation of this loop. 

We have calculated the mutation-induced stability changes 
of the ground-state conformations along the catalytic cycle of 
DHFR for the six distal mutations (see table|I|. These ground- 
state conformations correspond to the crystal structures of 
the five major intermediates for wildtype E. coli DHFR. We 
have used the program Concoord/PBSA ll60l . which is one 
of the most reliable programs for the calculation of mutation- 
induced stability changes I6T1I62") , and particularly suited for 
insertions of larger amino acids as in the mutation G121V 
since the flexibility of the protein backbone is taken into ac- 
count. A central assumption of our calculations is that the 
distal mutations only affect the conformational free energy 
of the enzyme but not the binding free energy of the ligands, 
which are "removed" from the crystal structures in the Con- 
coord/PBSA calculations. For the mutation G121 V, we obtain 
a destabilization of the closed ground-state conformations of 

and a stabilization of the oc- 



the intermediates E NH and E NH 



DHF' 



eluded ground-state conformations of the intermediates E^gp, 
Ejhf, and Ej^p with respect to the denatured state of the en- 
zyme (see table [I]). In agreement with the NMR experiments 
ll46l . our calculations thus indicate a stabilization of occluded 



conformations relative to closed conformations for G121V. 
These findings make qualitative sense in light of the protein 
structure. The larger Val side chain may interfere with the 
H-bond interactions between Asp 122 and Glyl5/Glul7 that 
would act to destabilize the closed conformation. For the mu- 
tation S 148 A, our calculations indicate a destabilization of the 
occluded ground-state conformations, and a rather neutral ef- 
fect on the closed ground-state conformations. The loss of the 
H bond between Serl48 and Asn23 would act to destabilize 
the occluded state, which is observed in the NMR experiments 

ma. 
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FIG. 7: Analysis of experimentally measured product unbinding 
rates k' o{f and calculated changes AAG co in conformational free- 
energy differences for six distal mutations of E. coli DHFR (see table 
2). The data are consistent with the conformational-selection route 
1 (red line with slope 1 according to eq. {TT}), but are not consistent 
with the induced-change route 2 (blue line, see eq. l|13[>). The data are 
also consistent with route 3 if the intermediate conformation of this 
route is similar to the closed conformation (see eq. |T7|). The data 
point for the mutant S148A (open circle) has to be excluded from 
this analysis since the closed conformation of this mutant seems to 
be the ground-state conformation prior and after product unbinding. 
In S148A, the product thus presumably unbinds directly from the 
closed conformation. 
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The mutation-induced change AAG co of the free-energy 
difference between the closed and occluded conformation can 
be calculated as the difference between the stability change of 
the closed ground-state conformation after product unbinding 
and the stability change of the occluded ground-state confor- 
mation in the product-bound state (see table 2). Within the 
calculation errors, the AAG co values are proportional to the 
values of In (k' oS / fc ff) obtained from the experimentally mea- 
sured off-rates of the wildtype and mutants (see fig. [7]), in 
agreement with eq. (JTTJ for the conformational-selection un- 
binding route 1, or with eq. ( 17 1 for route 3 and an intermedi- 



ate conformation i that is similar to the closed conformation. 
Our mutational analysis thus indicates that the conformational 
change from the occluded to the closed conformation occurs, 
at least to a large extent, prior to product unbinding. The mu- 
tation S148A (open circle in fig. |7|i has to be excluded from 
this analysis because the product presumably unbinds directly 
from the closed ground-state conformation of this mutant (see 
above). The mutation G121 V is included in our analysis since 
the effective unbinding rate along its conformational-selection 
route (see fig. [8]) still follows eq. |7} for route 1 if the relax- 
ation rate u' co for the additional conformational change into 
the occluded ground state after unbinding is much larger than 
the rebinding rate c' + [THF] in the closed conformation, which 
seems plausible. 



TABLE II: Mutational analysis of the DHFR product unbinding ki- 
netics 



mutation 


k'off 




-RT\n{k' oSS /k aS ) 


AAG co 


D122A 


3.5 ±0.1 


(i) 


0.75 ±0.10 


1.25 


D122N 


5.7 ±0.1 


(0 


0.46 ±0.10 


0.72 


D122S 


4.9 ±0.1 


« 


0.55 ±010 


0.27 


G121V 


1.9 ±0.3 


0) 


1.11 ±0.14 


1.98 


M42W 


175 




-1.56±0.10 


-1.42 


S148A 


113 ± 2 





-1.30 ±0.10 


-0.46 



w from ref. (57); (j) from ref. (49); (fe) from ref. (59|; (0 from 
ref. 1581 . - The values for — RT]n(k' g/k s) have been calcu- 
lated for the product unbinding rate k g = 12.5 ± 2 s -1 of 
the wildtype |48] and the given rates k' of{ of the mutants mea- 
sured at the temperature T — 25° G and pH=7. The mutation- 
induced changes of the free-energy difference between the closed 
and occluded conformation here is calculated from table 1 as 
AAG co = AAGco - AAGod where AAG oD = AAG(E NH ) 
is the mutation-induced stability change (free energy change rel- 
ative to the denatured state D) of the closed ground-state confor- 
mation after product unbinding, and AAG d = AAG(E^ H I F ) 
the stability change of the occluded ground-state conformation 
with bound product. The units of the k' o[f values are s~ , and 
the units of the values for — RT \n(k' of{ / k a f[) and AAG co are 
kcal/mol. The estimation of the error of AAG co is complicated 
by the fact that the mutation-induced free-energy change of the 
denatured state, which enters the AAG calculation, drops out 
in the difference since the denatured state is the same for both 
conformations. We assume that the error for AAG co is similar 
to the error for AAG and, thus, about 1 kcal/mol (see table 1). 
This assumption implies similar errors for the mutation-induced 
free-energy changes of the denatured state and the native confor- 
mations in the AAG calculations. 
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FIG. 8: Possible unbinding routes of the product THF from the 
E. coli DHFR mutant G121V, which has the occluded conforma- 
tion as ground-state conformation prior and after unbinding 1461 . 
The product thus may unbind from this mutant without conforma- 
tional change, or after a conformational change if the rate o'_ for 
direct unbinding from the occluded conformation is small. The 
conformational-selection route includes a change from the occluded 
to the closed conformation prior to product unbinding as on route 1 
of the wildtype (see fig.|4]l, followed by an additional conformational 
relaxation into the occluded conformation after unbinding. However, 
in analogy to eq. \\6) for route 3 of the wildtype, the effective un- 
binding rate along the conformational-selection route still follows 
eq. {7} for the wildtype route 1 if the relaxation rate u' co into the 
unbound, occluded ground state is much larger than the rebinding 
rate c' + [THF] in the closed conformation, which seems plausible. 



The extended catalytic cycle of wild-type DHFR 

NMR relaxation experiments indicate that wildtype E. coli 
DHFR populates a product-bound excited-state conformation 
that is rather similar to, but not identical with, the closed 
ground-state conformation after product unbinding [4 |. In ad- 
dition, the experiments indicate that the excited-state confor- 
mations in the product-bound state Ej^ F and in the unbound 
state E NH are similar [4], as on route (3) for product unbind- 



ing (see eq. ( 14 1). This structural information from the NMR 
relaxation experiments and our mutational analysis of the pre- 
vious section point towards a product unbinding route along 
which most, but not all, of the conformational change from 
the occluded to the closed conformation occurs prior to prod- 
uct unbinding. 

We focus here on the main conformational change from the 
bound ground state oEj^ F to the bound excited state cEj^f 
with a conformation similar to the closed conformation, and 
decompose the product unbinding process into two rather than 
three steps: 



6on 



£F NH 



cE^ 



(18) 



bgo c+[THF] 

The second step from the bound excited state cEj^ F to the un- 
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bound state cE NH includes a small conformational change be- 
sides the unbinding of the product. We do not further decom- 
pose this step since the conformational change during the step 
appears to be relatively small, and since the conformational 
transition rates in the unbound state have been determined in 
NMR relaxation experiments at the rather low temperature of 
8°C 0. The rates for the the main conformational transition 
in the product-bound state of the process ( 18 1 have been de- 



termined in NMR relaxation experiments at 27°C. These rates 
are b o£ = 18.5 ± 1.5 s" 1 and b £o = 735 ± 45 s" 1 0). 

The rates c_ and c + for the second step of the unbinding 
process ( 18 1 can be determined from the experimentally mea- 
sured values for the overall rates k s = 12.5 ± 2 s _1 and 
fe on = 2 ± 0.2 julVrV 1 at 25°C [48 1. The unbinding pro- 



cess ( 18 1 is similar to the conformational-selection unbinding 
route 1 (see eq. |6]i). From eq. |7]i for the effective off-rate of 
this route, we obtain 



b- co k s = 1500± 800 s -i 



k 



(19) 



off 



with the experimental values for and fc ot f given above. 

From eq. d8), we obtain 



{be 



ben 



6±2 /zM _1 s 



(20) 



with the experimental values for b C o and fc on and the value 
for c_ from eq. ( |19| l. The resulting value of c + justifies the 
assumption b^ 3> c + [THF] made in the derivation of eq. (JsJ 
since typical product concentrations [THF] are around 10 fiM 
in E. coli |63"1 . 

Our decomposition of the product unbinding step into the 
two substeps of eq . (p~8] > is included in the extended catalytic 
cycle shown in fig. 13] Together with the decomposition of the 
catalytic step in eq. the extended cycle specifies the or- 
dering of events in the coupling of the main conformational 
transition between the closed and occluded active-site loop 
conformations with binding and catalysis. 



DISCUSSION 

Role of the conformational changes in wildtype E. coli DHFR 

In contrast to E. coli DHFR, human and other vertebrate 
DHFRs only populate the closed conformation of the active- 
site loop along their catalytic cycle BTl l64l . This raises the 
question why there are conformational changes between a 
closed and an occluded active-site loop along the catalytic 
cycle of E. coli DHFR. If the conformational changes were 
'unfavorable', they could be 'eliminated' by mutations such 
as S148A that strongly destabilize the occluded conformation 
relative to the closed conformation. Similar to human DHFR, 
the mutant S148A of E. coli DHFR only populates the closed 
conformation (47), which is the required conformation for the 
chemical step. 

In contrast to human and other vertebrate cells, the internal 
pH of E. coli bacteria can vary from values below 5 ll65ll to 



peak values possibly larger than 8.7 [66], depending on the 
external pH of the surrounding medium. At low pH values, 
the overall catalytic rate of the extended catalytic cycle shown 
in fig.[3]is clearly limited by the conformational change prior 
to product unbinding. The conformational change with excita- 
tion rate 18.5 s _1 reduces the effective product unbinding rate 
to 12.5 s _1 (see eq. which is almost two orders of mag- 
nitude smaller than the maximal forward rate 950 s _1 of the 
chemical step attained at low pH values. At pH values close 
to the preferred internal pH 7.8 of E. coli ll65l . in contrast, 
the forward rate kf of the chemical step is of the same order 
of magnitude as the product unbinding rate (see eq. (|2]l). The 
change from the occluded to the closed conformation prior to 
product unbinding thus balances the strong pH dependence of 
the chemical reaction by limiting the overall catalytic rate at 
low pH to values of the same order of magnitude as the cat- 
alytic rate at the preferred pH of E. coli. 

Another difference between E. coli and human and 
other vertebrate cells concerns the relative concentrations of 
NADPH (NH) and NADP+ (N+). In prokaryotic cells such 
as E. coli, the concentrations of NH and N + are comparable 
||58l . In eukaryotic cells, in contrast, the concentration of N + 
is typically smaller than 1 % of the NH concentration. A con- 
sequence of the comparable concentrations of NH and N + is 
that E. coli DHFR works relatively close to the chemical equi- 
librium at pH values larger than the preferred pH value 7.8. In 
equilibrium, the product of all the rates in the forward direc- 
tion of the catalytic cycle is equal to the product of all rates in 
the backward direction ll67l . This equilibrium condition leads 
to 



1.3 • 10 n ~ pH [DHF] [NH] = [THF][N H 



(21) 



for the forward rate kf of the catalytic step given in eq. Q and 
the reverse rate k r = 0.55 s _1 at high pH. For the exemplary 
concentrations [DHF] = 0.3 /jM, [THF] = 13 fM, [NH] = 1 
mM, and [N + ] = 1.5 mM of Benkovic and coworkers ll63l . 



equilibrium is attained at pH = 9.3 according to eq. (21 1. The 
pH value at which equilibrium is attained is even lower for 
smaller substrate concentrations or higher product concentra- 
tion than in this example. Fluctuations of the substrate and 
product concentrations thus may lead to a reversal of the cat- 
alytic cycle of DHFR at pH values larger than the preferred 
pH of E. coli. 

In reverse direction, the catalytic rate of E. coli DHFR is 
limited by the low rate k r of the reverse catalytic step. On the 
extended catalytic cycle shown in fig. [3] the low value of k r 
is a consequence of the conformational change that precedes 
the chemical step in reverse direction. According to eq. Q, 
the actual rate of the reverse chemical step in the closed con- 
formation is fc_ ~ 38 s _1 at high pH, and the overall low re- 
verse rate k r ~ 0.55 s _1 results from the change between the 
occluded ground-state conformation and closed excited-state 
conformation with bound products THF and N + prior to the 
reverse chemical step. This conformational change thus limits 
the negative effects of a possible reversal of the catalytic cycle 
by strongly reducing the reverse catalytic rates. 

The conformational changes of E. coli DHFR thus may pro- 
vide robustness of the catalytic process against pH variations 
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and changes in the substrate, product and cofactor concentra- 
tions. In the forward direction of the catalytic cycle, the con- 
formational change prior to product unbinding prevents strong 
deviations of the overall catalytic rate for pH values below the 
preferred pH of E. coli. And at high pH, the conformational 
change during the catalytic step prevents high catalytic flux in 
the reverse direction. 



Stabilization of the occluded conformation can lead to low rates 
of the catalytic step 

Active-site remote mutations such as G121 V in the FG loop 
of E. coli DHFR (see fig. [TJ indicate that even distal residues 
in an enzyme can play important catalytic roles. The forward 
rate for the catalytic step of the mutant G121V is fcy ~ 1.4 s _1 
at pH 7 and, thus, significantly smaller than the wildtype 
rate kf ~ 220 s _1 at the same pH. We argue here that this 
strong decrease in the forward rate for the catalytic step is re- 
lated to a stabilization of the occluded conformation relative to 
the closed conformation. NMR experiments with the mutant 
G121 V indicate that the occluded conformation is the ground- 
state conformation in all five intermediates along the catalytic 
cycle [46 1 . Since hydride transfer from NH to DHF is only 
possible in the closed conformation, the catalytic step of the 
mutant G121V requires an excitation into the closed confor- 
mation prior to the chemical reaction, and a relaxation into the 
occluded conformation after the reaction: 



the mutant G121V is the product of the equilibrium constant 
K' s and the rate k' + for the chemical substep. In contrast, the 
forward rate k t for the catalytic step of the wildtype is identi- 
cal with the rate kf of the chemical substep (see eq. |5])). We 
would like to emphasize that eq. ( 23 1 is in agreement with the 
observed kinetic isotope effect for the forward rate of G121 V 
ll49l |6"8if70l since the rate k'f is proportional to the rate k' + 



of the actual chemical step. The catalytic scheme (22 1 is also 
in agreement with pre-steady-state kinetic experiments, which 
indicate a conformational change prior to the chemical reac- 
tion with an excitation rate of s' oc ~ 3.5 s _1 11491 . Similar 
to G121V, the strong reduction in the forward rate of the cat- 
alytic step observed for other distal mutations such as D122A, 
D122N, and D122S |[57l may also results from a stabilization 
of the occluded conformation. 

We have attributed here strong reductions in the forward 
rate of the catalytic step caused by distal mutations to changes 
in the conformational equilibrium. For mutations within the 
active-site loop (Met20 loop), a different, dynamic mecha- 
nism has been recently proposed by Bhabha et al. J47). The 
mutant N23PP and the double mutant N23PP/S148A only 
populate the closed state of the active-site loop, but reduce 
the forward rate kt for the catalytic step by a factor of 15. 
This reduction is traced back to a stiffening of the active- 
site loop by the mutation N23PP, which dampens loop fluc- 
tuations that seem to be required for catalysis in the closed 
conformationflTl . 



oE 



NH 
DHF 



cE 



NH 
DHF 



k' 



N+ 
THF 



1J 



oE: 



N+ 
THF 



(22) 



CONCLUSIONS 



In contrast, the catalytic step of the wildtype only involves 
a conformational relaxation after the chemical substep (see 
eq. Q). The effective forward rate for this catalytic 3-step 



process of the mutant G121V is (see eq. ( 39 1) 



k'r 



for k' > k'_ 



^k' + = K' s k' + for4,»fcV 



(23) 



with the constant K' s = s' oc /s' co for the conformational equi- 
librium of the substrate-bound states. We have assumed here 
that the relaxation rate k' co into the product-bound ground 
state oEj^j F is much larger than k'_, as for the wildtype (see 
fig. 13}. In addition, we have assumed that the relaxation rate 
s' co into the substrate-bound ground state is much larger than 
k',, which is reasonable at least for intermediate pH values 
(see eqs. |2} and {5}). 

The small forward rate fcy for the catalytic step of the mu- 
tant G121V is a consequence of the fact that the conforma- 
tional equilibrium constant K' s in eq. ( 23 i is much smaller 



than 1 since the occluded conformation is the ground-state 
conformation. According to eq. ( 23 1, the forward rate k'f of 



In this article, we have presented extended catalytic cycles 
for the enzyme DHFR from E. coli that specify the order- 
ing of chemical and physical steps, including conformational 
changes. The extended cycles were derived in a general the- 
oretical framework that provides a link between the confor- 
mational transition rates from dynamic NMR experiments and 
the effective binding and catalytic rates measured in relaxation 
experiments. Our theoretical framework allows to distinguish 
between induced conformational changes that occur predomi- 
nantly after a binding/unbinding or chemical process, and se- 
lected conformational changes that occur predominantly prior 
to such processes. In addition, the framework helps to explain 
the effect of mutations on the rates of binding and catalysis. 
Our approach is rather general and applicable to other pro- 
teins and experimental methods that probe the conformational 
or binding dynamics. 



APPENDIX 

Relaxation rates of the three-state process 

The catalytic process of eq. Q and the routes 1 and 2 for 
product unbinding of E. coli DHFR follow the general 3-state 
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reaction scheme 



with P(t) = (P A (t),P B (t),P c (t)) and 



kAB k BC 

A^± B ^± C 

&BA k C B 



(24) 



The time-dependent probabilities P A (t), P B (t), and Pc{t) of 
the three states A, B, and C in this scheme are governed by 
the master equations 



dP A (t) 

dt 
dP B (t) 

dt 
dPcjt) 

dt 



k BA P B (t) - k AB P A (t) 

k AB P A {t) - (k BA + k BC )P B (t) + k CB P c {t) 
k BC P B {t) - k CB P c (t) 



where denotes the rate from state i to state j. The three 
equations can be written in the matrix form 



dP(t) 
dt 



= -WP(t) 



(25) 



W 



k AB ~k BA 
-k AB k BA + k BC -k CB 
-k BC k CB 



The general solution of eq. (25 1 has the form [71 ] 



P(t) = c Y + c 1 Y 1 e 



c 2 Y 1 e 



-k 2 t 



(26) 



(27) 



where Y a , Y\, and Y2 are the three eigenvectors of the ma- 
trix W, and ki and k 2 are the two positive eigenvalues 



k\ = \ (k AB + k BA + k BC + k CB - \J {k AB + k BA + k BC + k CB ) 2 - 4:(k AB (k BC + k CB ) + k CB k BA ) 



ko = 



^ (kAB + k BA + k BC + k CB + \J (k AB + k BA + k BC + k CB ) 2 - 4:(k AB (k BC + k CB ) + k CB k BA )^ 



(28) 



(29) 



of W with corresponding eigenvectors Y\ and Y 2 - The 
smaller eigenvalue k\ is the dominant relaxation rate on long 
time scales. This eigenvalue corresponds to the relaxation rate 



ka 



fei 



(30) 



observed in experiments. 



Effective rates of the three-state process 

A general effective rate k A ^c from state A to state C can be 
defined by considering an adsorbing state C without probabil- 
ity backfiow into state B. An absorbing state C with kc B = 
corresponds, e.g., to a situation in which the probability flow 
from state C into another state on a chemical cycle is much 
faster than the backfiow into B. The effective rate from A to C 
is the dominant relaxation rate k\ for kc B =0: 



kA-yC = h 



k C B=0 



(31) 



The effective rate k A ^c depends on the three rates k AB , k BA , 



and k B c (see eq. (28 1). In limiting cases in which one of these 



rates is much smaller than another rate, we obtain 

k AB k BC 



k A - 



k BA + k BC 
k AB k BC 



for k AB < k BA or k AB < k B c 
for k BC < k BA or k BC < k AB 



k AB + k BA 

m.in(k AB , k BC ) for k BA < k AB or k BA < k BC 

(32) 

An important special case of the process ( |24| is the case 
in which the equilibrium probability of state B is significantly 
smaller than the equilibrium probabilities of A and C. In this 
case, we have k AB <C k BA and kcs <S k B c, and the relax- 
ation rate <f30b is 



(33) 



k AB k BC + k BA k CB . , , 

Kobs — ; — ; = k A ^ c + kc^rA 

k B A + k BC 

with the effective rates 

kAsk B c 



k A - 



■c 



and 



^C-¥A 



k B A + k B c 



k B AkcB 
k B A + k B c 



for k AB < k BA 



for k CB < k BC 



(34) 



(35) 



obtained from eq. ( 32 1. In the special case of a low -probability 
intermediate state B, the dominant relaxation rate d33l is also 



12 



obtained in the steady-state approximation for the 3-state pro- 
cess ( |24| . The steady-state approximation is based on the as- 
sumption dP B (t)/dt = 0, i.e. on a negligible change of the 



probability for the state B in the set of equations (25 i to (25 1. 
In general, however, the dominant relaxation rate k a \, s is not 
equal to k A ^c + k c ^A- 

Our derivation of the eqs. Q, Q, Q, and ^ for the effec- 
tive rates of the catalytic step and the product unbinding step 



along the catalytic cycle of E. coli is based on eq. ( 32 1 and, 
thus, on relaxation rates. This derivation is in line with the 
enzyme kinetics measurements, in which the rates for these 
steps are determined from pre-steady-state or relaxation ex- 
periments ll48l . Alternatively, the effective rates for these 
steps can be derived from results for the steady-state rate of 
the enzymatic cycle (T. R. Weikl, unpublished). 



Effective rates of a four-state process 

Route 3 for product unbinding of E. coli DHFR follows the 
general four-state reaction scheme 



kAB k BC k C D 

A ^± B ^± C ^± D 

kBA kcB kr>c 



(36) 



easily, which leads to the relaxation rate 

, ^ k AB k BC k C D + k BA k CD k DC 
k B AkcB + k B AkcD + ksckcD 



(38) 



For kAB <C k B A and kuc "C &cd, the approximate expres- 



sion (38 1 is an excellent agreement with numerical results for 



the dominant relaxation rate of the process ( 36 1, which is the 
smallest nonzero eigenvalue of the matrix ( 37 1. 

Similar to the three-state process with low -probability inter- 
mediate B, the relaxation rate fc bs is the sum of the effective 
rate 



k A - 



kABkBckcD 



ksAkcB + k B AkcD 
from A to D and the effective backward rate 

k B AkcokDC 



kncX 



BC^CD 



'"D^A 



ksAkcB + ksAkcD + k B ckcD 



(39) 



(40) 



from D to A. The effective rate form A to D is the dominant 
relaxation rate for kcD — 0, i.e. for an absorbing state D 
without probability backfiow into C. The effective rate from 
D to A is the dominant relaxation rate for k AB = 0, i.e. for an 
absorbing state A. 



The probability evolution of the four states is governed by the 
master equation set 

dP A (t) 

dt 
dP B (t) 

dt 
dPcjt) 

dt 
dP D {t) 

dt 

These equations can be written in the form of the matrix equa- 
tion d25l with P(t) = {P A {t),P B (t),P c (t),P D {t)) and 



= k BA PB{t)~ 


k AB PA{t) 




= k AB P A {t)- 


(k B A + k BC )PB(t)^ 


- k CB Pc(t) 


= k B cP B {t)- 


{k C B + k CD )Pc(t)- 


^k DC P D {t) 


= k CD Pc(t)- 


k D cP D (t) 





w 



( k AB -k BA \ 

-k AB k BA + k B c -kcB 

-fcfic k C B + k C D -k DC 

V -k CD k DC ) 



(37) 



The three relaxation rates of the four-state process are the 
three positive eigenvalues of the matrix W. 



The states B and C of the process ( 36 1 correspond to the 
excited states iEj^ F and iE NH along route 3, which exchange 
with the ground states with rates b Q i <C bi and u c j <C Ui C . 
We focus here on the corresponding case k A s <C ks A and 
koc "C kcD- m this case, the intermediate state B has a 
low equilibrium probability compared to A, and the state C a 
low probability compared to D. This makes the steady-state 
approximation plausible, which assumes dPs(t)/dt = and 
dPc(t)/dt — for the intermediate states. Within this ap- 



proximation, the four equations ( 37 1 to ( 37 1 can be solved 
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